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Abstract 

This paper presents the experimental and 
theoretical results for expanding methane and ethylene 
diffusion flames in microgravity. A small porous sphere 
made from a low-density and low-heat-capacity insulating 
material was used to uniformly supply fuel at a constant 
rate to the expanding diffusion flame. A theoretical 
model which includes soot and gas radiation is formulated 
but only the problem pertaining to the transient expansion 
of the flame is solved by assuming constant pressure 
infinitely fast one-step ideal gas reaction and unity Lewis 
number. This is a first step toward quantifying the effect 
of soot and gas radiation on these flames. The 
theoretically calculated expansion rate is in good 
agreement with the experimental results. Both 
experimental and theoretical results show that as the flame 
radius increases, the flame expansion process becomes 
diffusion controlled and the flame radius grows as Vl 
T heoretical calculations also show that for a constant fuel 
mass injection rate a quasi-steady state is developed in the 
region surrounded by the flame and the mass flow rate at 
any location inside this region equals the mass injection 
rate. 

I. Introduction 

The absence of buoyancy-induced flows in a 
microgravity environment and the resulting increase in the 
reactant residence time significantly alters the 
fundamentals of many combustion processes. Substantial 
differences between normal gravity and microgravity 
flames have been reported during droplet combustion! 1], 
flame spread over solids[2,3], candle flames[4] and others. 
These differences are more basic than just in the visible 
flame shape. Longer residence time and higher 
concentration of combustion products create a 


thermochemical environment which changes the flame 
chemistry. Processes such as soot formation and 
oxidation and ensuing flame radiation, which are often 
ignored under normal gravity, become very important and 
sometimes controlling. As an example, consider the 
droplet burning problem. The visible flame shape is 
spherical under microgravity versus a teardrop shape 
under normal gravity. Since most models of droplet 
combustion utilize spherical symmetry, excellent 
agreement with experiments is anticipated. However, 
microgravity experiments show that a soot shell is formed 
between the flame and the evaporating droplet of a sooty 
fuelfl]. This soot shell alters the heat and mass transfer 
between the droplet and its flame resulting in significant 
changes in the burning rate and the propensity for flame 
extinction. This change in the nature of the process 
seems to have occurred because of two reasons: (i) The 
soot formed could not be swept out of the flame due to 
the absence of buoyant flows. Instead, it was forced to go 
through the high temperature reaction zone increasing the 
radiative heat losses, and (ii) soot formation was enhanced 
due to an increase in the reactant residence time. 

Recently, some very interesting observations of 
candle flames under various atmospheres in microgravity 
have been reported[4]. It was found that for the same 
atmosphere, the burning rate per unit wick surface area 
and the flame temperature were considerably reduced in 
microgravity as compared with normal gravity. Also, the 
flame (spherical in microgravity) was much thicker and 
further removed from the wick. It thus appears that the 
flame becomes "weaker” in microgravity due to the 
absence of buoyancy generated flow which serves to 
transport the oxidizer to the combustion zone and remove 
the hot combustion products from it. The buoyant flow, 
which may be characterized by the strain rate, assists the 
diffusion process to execute these essential functions for 
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the survival of the flame. Thus, the diffusion flame is 
"weak” at very low strain rates and as the strain rate 
increases the flame is initially "strengthened" and 
eventually it may be "blown-out." The computed 
flammability boundaries[5] show that such a reversal in 
material flammability occurs at strain rates around 5 sec' 1 . 
Model calculations for a zero strain rate 7-D diffusion 
flame show that even gas radiation is sufficient to 
extinguish the flame [6]. 

The above observations suggest that flame 
radiation will substantially influence diffusion flames 
under microgravity conditions, particularly the conditions 
at extinction. This is because, flame radiation at very low 
or zero strain rates is enhanced due to: (i) high 
concentration of combustion products in the flame zone 
which increases the gas radiation, and (ii) low strain rates 
provide sufficient residence time for substantial amounts 
of soot to form which is usually responsible for most of 
the radiative heat loss. It is anticipated that this radiative 
heat loss may extinguish the already "week" diffusion 
flame. 

To investigate the possibility of radiation-induced 
extinction limits under microgravity conditions, spherical 
geometry is chosen. This is convenient for both 
experiments and theoretical modeling. In this work, a 
porous spherical burner is used to produce spherical 
diffusion flames in pg. Experiments conducted with this 
burner on methane (less sooty) and ethylene (sooty) 
diffusion flames are described in the next section. A 
general theoretical model for transient radiative diffusion 
flames is then formulated and calculations are presented 
for the transient expansion of the spherical diffusion 
flame. These calculations are compared with the 
experimental measurements in the discussion section. 
This work is the first necessary step toward investigating 
radiative-extinction of spherical diffusion flames. 

II. Experimental Apparatus and Results 

The pg experiments were conducted in the 2.2 
sec drop tower at the NAS A Lewis Research Center. The 
experimental drop-rig used is schematically shown in 
Color Plate 1. It consists of a test chamber, burner, 
igniter, gas cylinder, solenoid valve, camera, computer 
and batteries to power the computer and the solenoid 
valves. The spherical burner (1.9 cm in diameter) is 
constructed from a low density and low heat capacity 
porous ceramic material. A 150 cc gas cylinder at 
approximately 46.5 psig is used to supply the fuel to the 
porous spherical burner. Typical gas flow rates used were 
in the range of 3-15 cm 3 /s. Flow rates to the burner are 
controlled by a needle valve and a gas solenoid valve is 
used to open and close the gas line to the burner upon 


computer command. An igniter is used to establish a 
diffusion flame. After ignition the igniter is quickly 
retracted from the burner and secured in a catching 
mechanism by a computer-controlled rotary solenoid. 
This was necessary for two reasons (i) The igniter 
provides a heat sink and will quench the flame (ii) Upon 
impact with the ground (after 2.2 sec) the vibrating igniter 
may damage the porous burner. 

As shown in the Color Plate 1, the test chamber 
has a 5" diameter Lexan window which enables the 
camera to photograph the spherical diffusion flame. The 
flame growth can be recorded either by a 16mm color 
movie camera or by a color CCD camera which is 
connected to a video recorder by a fiber-optic cable 
during the drop. Since the fuel flow may change with 
time, it had to be calibrated for various settings of the 
needle valve for both methane and ethylene. A soap 
bubble flow meter was used to calibrate the flow for 
various constant gas cylinder pressures. Constant 
pressures were obtained by connecting the cylinder to the 
main 200 lb gas cylinder using a quick-disconnect An 
in-line pressure transducer was used to obtain the transient 
flow rates. Changes in the cylinder pressure during the 
experiment along with the pressure-flow rate calibration, 
provides the transient volumetric flow rates. These are 
shown plotted in Figure 1. 



Figure 1: Volume flow rales versus time. 


In Figure 1, the letters "M" and "E" represent 
methane and ethylene respectively and the letters "L", 
"M" and "H" represent low, medium and high flow rates. 
Thus, MM implies medium flow rate of methane. Note 
that low flow rate for methane is nearly equal to the 
medium flow rate of ethylene. For these experiments, the 
gas velocity at the burner wall was between 0.25-1 
cm/sec. 

The porous spherical burner produced a nearly 
spherical diffusion flame in microgravity. Some observed 
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Picture of the microgravity spherical diffusion flame apparatus 
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disturbances are attributed to slow large-scale air motion 
inside the test chamber. Several microgravity experiments 
were performed under ambient pressure and oxygen 
concentration conditions for different flow rates of 
methane and ethylene (as shown in Fig. 1). Methane was 
chosen to represent a non-sooty fuel and ethylene was 
chosen to represent a moderately sooty fuel. In these 
experiments, ignition was always initiated in 1-g just prior 
to the drop. The package was typically dropped within 
one second after ignition. The primary reason for not 
igniting in pg was the loss of time in heating the igniter 
wire and in stabilizing the flame after the initial ignition 
disturbances. Some photographs from these experiments 
are shown in the Color Plate 2. 

The flame radius measured from such 
photographs along with the model predictions (to be 
discussed later) are shown in Fig. 2. As expected, for the 
same flow rates it was found that ethylene flames were 
much sootier and smaller. Immediately after dropping the 
package, the flame shape changed from a teardrop shape 
(see Color Plate 2) to a spherical shape (although it was 
not always completely spherical, probably because of slow 
large-scale air motion persisting inside the test chamber). 
The photographs shown in the Color Plate 2 are for 
medium flow rates of methane and low flow rates of 
ethylene. For the data presented in Fig. 2, an average 
flame radius determined from the photographs was used. 
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Figure 2: Flame radius versus time 


It is interesting to note that for both methane and 
ethylene (see the progressive flame growth in Color Plate 
2), initially and in 1-g (e.g. photographs ’e’ & T) the 
flame is nearly blue (non-sooty) but becomes bright 
yellow (sooty) immediately after the onset of pg 
conditions. Later, as the pg time progresses, the flame 
grows in size and becomes orange and less luminous and 
the soot seems to disappear. A possible explanation for 
this observed behavior is suggested by the theoretical 
calculations of Ref. 6. The soot volume fraction First 
quickly increases and later decreases as the local 
concentration of combustion products increases. 
Essentially, further soot formation is inhibited by the 
increase in the local concentration of the combustion 
products [Ref.7,8] and soot oxidation is enhanced. Thus, 
at the onset of pg conditions, initially a lot of soot is 
formed in the vicinity of the flame front (the outer faint 
blue envelope) resulting in bright yellow emission. As 
the flame grows, several events reduce the flame 
luminosity: (i) The soot is pushed toward cooler regions 
by thermophoresis. In fact, for sootier fuels this leads to 
the formation of a soot shell, (ii) The high concentration 
of combustion products left behind by the flame front 
inhibits soot formation and promotes soot oxidation, (iii) 
The dilution and radiative heat losses caused by the 
increase in the concentration of combustion products 
reduces the flame temperature which in turn reduces the 
soot formation rate and the flame luminosity. 

Figure 2 shows the average measured flame 
radius for methane and ethylene pg diffusion flames 
plotted against time. This is the radius of the outer faint 
blue region of the flame as measured from the 
photographs. To a good approximation this may be 
considered as the flame front location. Thus, as a first 
step, it will be interesting and important to determine if 
the transient expansion of the pg spherical diffusion flame 
can be theoretically predicted without considering soot 
formation and oxidation kinetics and flame radiation. 

III. Model Formulation 


As noted above, the spherical diffusion flames 
are expanding and changing their luminosity with time. 
Thus, the general theoretical formulation must be transient 
and must include flame radiation. For the simplest case 
of constant pressure ideal gas reactions with Le=l, we 
may write the following governing equations for any 
geometrical configuration (spherical or counterflow 
geometry): 

Mass Conservation: 


.^.+V<pv)=0 


( 1 ) 
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Energy Conservation : 


Oxygen Conservation: 


p^L+pv-Vh ‘-V-ipDVh 0 
ot 

=-Z h :w r Qmr-v-<r, 

i 

Constant Pressure Ideal Gas: 

pT=pX or ph ' = const. 


( 2 ) 


(3) 


Here, the symbols have their usual definitions 
with p = density, T = temperature, v = velocity, = 
mass fraction of species i, h s = sensible enthalpy, Wj = 
mass production or destruction rate per unit volume of 
species i and D = diffusion coefficient The last three 
terms in Equ (2) respectively are: the chemical heat 
release rate due to gas phase combustion, chemical heat 
released due to soot oxidation and the radiative heat loss 
rate per unit volume. The above equations, however, are 
insufficient for our problem because the soot volume 
fraction must be known as a function of space and time 
to determine the radiative heat loss. To enable describing 
soot volume fraction in a simple manner, we define the 
mass fraction of atomic constituents as follows: 




, where M* is the molecular weight 


of species i, is the atomic weight of atom j and v/ is 
the number of atoms of kind j in specie i. Assuming that 
the only atomic constituents present in the hydrocarbon 
flame are C, H, O & Inert and with = d> = p s f v /p 
(where: p s = soot density & f v = soot volume fraction), 


we obtain: C c + * pjjp = 1 • 

Defining 4c + 4 h = 4f and Zp = 4 f /Y f „ and Z 0 = 4o 


we obtain Z=[($ F ) F „Z F + pjjp] as the 

conserved scalar for a sooty flame. This yields the 
following soot, fuel and oxidizer conservation equations 
in terms of their scalar variables: 


Soot Conservation: 
54> 


p +p v-V(<!>) -V ip DV(<P)] 
ot 


-mf-mf = m~, 


(4) 


Fuel Conservation: 


p^L + p^V(Z F )-VipDV(Z F )) 

ot 

- -—(m S ' m s )- — - 

1 Foo 1 F- 


(5) 


P 


dZ 0 


~5T 


+pPV(ZJ-V1pDV0J]-O 


( 6 ) 


Under conditions of small soot loading, the soot 
terms in the fuel and energy conservation equations can 
be ignored except when studying radiative extinction. 
Thus, Equ (5) may be considered homogeneous to a good 
approximation. Also, as a first crude approximation, the 
heat lost by flame radiation may be subtracted from the 
heat of combustion in the form of a radiative fraction. 
Thus, the energy equation (Equ(2)) can also be made 
homogeneous if written in terms of the total enthalpy [h 
= Zj Yj (l^ 0 + hi 5 )]. This approach may be adequate for 
calculating the observed expansion rate of the spherical 
diffusion flames, but it is completely inadequate for 
predicting radiative extinction. However, the great 
mathematical advantage of this approach is that it makes 
Eqs. (2, 4, 5, & 6) identical and only one conserved scalar 
equation need be considered. As a first step, it is of 
interest to see how well the transient expansion of the pg 
spherical diffusion flames be predicted without rigorously 
considering soot and gas radiation. This will also help in 
quantifying the effect of soot and gas radiation by 
comparison with more detailed calculations. Re-writing 
the above equations in spherical coordinates, we get: 


Mass Conservation: 




Fuel Conservation: 

dZ dZ 

p ir +pv V 


1 d f . n dz 

r-pD 

dr 


= 0 


(7) 


( 8 ) 


These two equations along with the ideal gas law 
at constant pressure, Equ.(3), are sufficient to describe the 
transient growth of non-radiative spherical diffusion 
flames and are expected to approximate this growth in the 
presence of flame radiation. It is also assumed that a fast 
one-step overall reaction occurs at the flame surface. This 

N-2 

is represented by: y F + V O — > v P \ with 

i»l 

q° as the standard heat of reaction and Q = q7M F v F the 
heat released per unit mass of fuel. Clearly, 

q » = H>e 

1*1 

corresponding initial and boundary conditions for a sphere 

of radius ’R’ blowing fuel gases at a rate A/(f)are 
discussed below and illustrated in Figure 3. 
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Figure 3: Schematic of the Model Problem 


Continuity, fuel mass fraction and energy 
conservation at the surface of the sphere yield: 


M(t) 
4nR 2 


=(p v L 


(9a) 


M(t) 

4kR 2 


(Y f .-Y„) = -PD 


f \ 

M F 

IF 


>* 


(9b) 


m_, hL 

4kR 1 


-V)=-pz> 


dh 

sr 


>/? 


(9c) 


Here, Y Foo & are the fuel mass fraction and enthalpy 
of the incoming fuel stream and Yp* and h s R are the 
corresponding values at the outer surface of the sphere. 
The ambient values of fuel and oxidizer enthalpies are 

taken to be equal i.e. /ij = and Z=1 & 

p^p 0 for ambient conditions on the fuel side and ZM) & 
p=p 0 for ambient conditions on the oxidizer side. Now, 
for high fuel injection rates, Yp- ~ Yfr ; h^ ~ h s R and 
Zr* 1 and the corresponding diffusion terms in Equs. 9b 
& 9c become zero. For a given mass injection rate [ 

M(t) ]> these conditions are also satisfied as R— >0. Thus, 
for a point fuel source, the boundary conditions at the 
source are simplified. Other initial and boundary 
conditions are: At t=0, Z(r,0); p(r,0) & v(r,0) are the 
spatial distributions corresponding to the flame at t=0, as 
shown in Figure 3. Also, at the flame surface [r^t)] 

Z-Z e .(l -Y F _M 0 v 0 /r o _M F v F r'. and aa r-*~. 


Z->0; v-*0 & h s — >h s „. All other variables can be easily 
obtained in terms of Z by utilizing the linear relationships 
between the conserved scalars [Ref. 10]. For constant 
pressure ideal gas reactions these linear relationships 
yield: 


For R < r < r f (t): 

p = Pc 




h:a-z) 


J 


( 10 ) 


m) _pd P 

4jtr 2 p o p Ft 


+ , 


M(t) QZ C 
4*r 2 p o hf 



y ( Z *~ Z ) 

f -(i-z; 


(ii) 


These equations, along with Equ. (8), are 
sufficient to provide all the distributions in the region 
between the porous sphere and the flame. In Equ. (11), 
the first term on the right hand side represents the 
injection velocity and the second term accounts for the 
increase in the velocity due to the decrease in density. 
The third term is identically zero if the distribution of Y F 
within the porous sphere (r<R) is identical to that in the 

gas i.e. Y f = Y f JZ-Z c )I (1 -Z c ) . Note at n=R, Y F 
= Y re and Z = Zr. Also note that the third term becomes 
zero for high injection velocities and small ’R’ since 
Zr-> 1 & YpR— >Y Foo . In Equ. (11), p and (dp / dr) can 
be expressed entirely in terms of Z through Equ. (10). 
Thus, Equ. (8) along with the appropriate boundary 
conditions is sufficient to determine Z(r,t). 


At the flame surface r=rft): 


At the flame surface, the Z ( ‘ } = Z ( *° = Z c and all 
its derivatives are continuous. Here, represents the 
fuel side and represents the air side. Also, T(r/\ t) = 
T(r f \ t) = T f ; p(r f + , t) = p(r f \ t) = p f and v(r f \ t) = v(r f \ 
t) = v f . Other jump conditions at the flame surface are 
obtained from species and energy balances as follows 
(assuming Le = 1 & D = D t ): 


f \ 


( \ 


M f 

= _ 

M f v f 

M 0 

l dr J 

r,J) 

l W ° V °J 



dh 

“5T 




r dh’^ 
IF 


= -Q 




r \ 

~FF 


Mr A 


03) 


In terms of Z, both Equs. (12) & (13) are 
identically satisfied if the first derivatives of Z are equal 
at the flame surface. Thus, for the solution of Equ. (8) in 
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the domain r > r f , we only need to find expressions for p 
and v in terms of Z. 


For r f (t) < r < oo; 

P = Pc 




hi 


j 


(14) 




m 

' (\~Z) 4 7tr/p 

M(t)QZ c 

^r/ph: 


f 

f z R -z Y 

Y-Y c 


FR F<* 

1-Z 

\ 

c )) 


]- Dd P (15) 
J P "IF 


In the derivation of Equ. (15), gas velocity and 
density at the flame front are made continuous i.e. v(r f \ 
t) = v(r f \ t) = v f . and p(r f +, t) = p(r f \ t) = p f . Thus, v f in 
Equ. (15) can be obtained from Equ. (11). Once again, 
the third term inside the bracket of Equ (15) becomes 
zero for reasons discussed above. Equs. (14), (15) and 
(8) along with the boundary conditions are sufficient to 
determine Z(r,t) for n>r f . 

IV Solution 


Before discussing the solution procedure, let us 
examine the porous sphere used in the experiments. This 
sphere is quite small (19 mm dia.) and is constructed 
from a high porosity, low density and low heat capacity 
insulating material. Thus, its capacity to store heat and 
mass is negligible compared to the fuel injection rate 
which is injected inside the sphere (see Fig. 3). Hence, 
conditions inside the material of the sphere equilibrate on 
a time scale much shorter than the flame expansion time 
i.e. convection balances diffusion for any variable under 
consideration that can be described by an equation similar 
to Equ. (8). Neglecting radiation from the surface of the 
sphere, conservation conditions yield equations identical 
to Equs. 9(a) and 9(b) where ambient conditions are 
assumed to exist near the center of the sphere. 
Physically, the only purpose the porous sphere serves is 
to provide a radially uniform flow and it does not 
participate in energy and species balances because of its 
low storage capacity. Thus, the boundary conditions at 
the source can be applied at an arbitrarily small radius ’R’ 
(chosen for numerical convenience) such that Y Foo = 
Yf*; h 5 ^ « h s R and Z* ~ 1. This considerably simplifies 
Eqs. (11) & (15). 


Equation (8), with Eqs. (10) & (11) for the fuel 
side and Eqs. (14) & (15) for the air side were 
numerically solved using the method of lines. A 
computer package entitled DSS2 was employed for this 
purpose. The calculated results for the flame location are 
shown plotted by dotted lines in Fig. 2. Property values 
used were those for air (p c =1.16 x 10' 3 gm/cm 3 , D 0 = 
0.226 cm 2 /s, X. =298K, Cp =1.35 J/kgK) and the 
diffusion coefficient was assumed to vary as T 3/2 as 
predicted by kinetic theory of gases. Heat of combustion 
(Q) and mass based stoichiometric coefficient (v) used 
for methane and ethylene were Q=47465 J/gm and v=4 
and Q=47465 J/gm and v=3.429 respectively. No 
assumptions other than those stated above were made to 
match the experimental data. Initial spatial distribution of 
Z(r,0) required for the flame at the start of pg time (i.e. at 
t=0) was taken as: 


Z(r, 0) = erfc 


0 r-R)erfc-'(Z t ) 

(P^R) 


(16) 


V Results and Discussion 


Figure 2 shows the average radius of the outer 
faint blue regions for both methane and ethylene pg 
diffusion flames plotted against time. This radius was 
measured from the photographs. As stated above, the 
corresponding calculated results for the flame location are 
shown plotted with dotted lines. Given the 

approximations made in the model and the experimental 
errors, the comparison between the experimental and 
predicted flame radius is quite encouraging. 

Numerical calculations also yield the 

instantaneous velocity and density profiles around the 
porous sphere during the flame expansion. These are 
shown plotted in Figures 4 & 5. Starting from the porous 



Figure 4: Radial velocity distribution at various instants 
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sphere (r=0.95), the gas velocity drops sharply and 
becomes a minimum at the flame location (r=r f ). 
Surprisingly, the mass flow rate at any location r < r f is 
found to be equal to the mass injection rate (i.e. 

47 T p vr 2 =M(t) )• implies that a similarity exists in 

the normalized coordinate r/r^t) in the region r < r^t). 
The density profiles in this region (Fig. 5) also show a 
similarity. Further reflection shows that this is to be 



Figure 5: Radial density distribution at various instants 

expected. In this problem, a constant temperature 
(adiabatic flame temperature) spherical flame is 
propagating outward starting from a small radius. In the 
spherical geometry, heat loss from the region surrounded 
by the flame is not possible. Thus, the only heat required 
by this region from the flame (in the absence of radiation) 

is to heat the injected mass M(f) 10 flame 

temperature. Since, the injected mass is taken to be 
constant with time, a quasi steady state is developed. 
This is also observed in the density gradients at the flame 
on the fuel side (which are constant and are proportional 
to the temperature gradients). Applying a simple energy 
balance over the region r < r f , we obtain: 

M(t) _ M(t) (1?) 

4jtr/p / 4nrf p o P, J. ' 

Using Equ. (11) we find that at the flame 

Vj. = M(t)/4n pj. ^ important to note that this 

is possible only bei^use the injection rate is not varying 
with time. 

VI Conclusions 

In this work, experimental and theoretical results 
for expanding methane and ethylene diffusion flames in 
microgravity are presented. A small porous sphere made 
from a low-density and low-heat-capacity insulating 
material was used to uniformly supply fuel at a constant 


rate to the expanding diffusion flame. A theoretical 
model which includes soot and gas radiation is formulated 
but only the problem pertaining to the transient expansion 
of the flame is solved by assuming constant pressure 
infinitely fast one-step ideal gas reaction and unity Lewis 
number. This is a first step toward quantifying the effect 
of soot and gas radiation on these flames. The 
theoretically calculated expansion rate is in good 
agreement with the experimental results. Both 
experimental and theoretical results show that as the flame 
radius increases, the flame expansion process becomes 
diffusion controlled and the flame radius grows as Vt. 
Theoretical calculations also show that for a constant fuel 
mass injection rate a quasi-steady state is developed in the 
region surrounded by the flame and the mass flow rate at 
any location inside this region equals the mass injection 
rate. 
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